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We study the properties of quantum single-particle wave pulses created by sharp-edged or apodized 
shutters with single or periodic openings. In particular, we examine the visibility of diffraction fringes 
depending on evolution time and temperature; the purity of the state depending on the opening-time 
window; the accuracy of a simplified description which uses "source" boundary conditions instead 
of solving an initial value problem; and the effects of apodization on the energy width. 
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I. INTRODUCTION 

Diffraction in time was discussed first by Moshinsky 
The hallmark of this phenomenon consists in tem- 
poral oscillations, deviating from the classical regime, 
of Schrodinger, matter waves released in one or several 
pulses from a preparation region in which they are ini- 
tially confined. The original setting consisted of a sud- 
den opening of a shutter to release a semi-infinite beam, 
and provided a quantum, temporal analogue of spatial 
Fresnel-diffraction by a sharp edge Q ; later more compi- 
cated shutter windows, initial confinements, and time-slit 
combinations have been considered , in particular two 
temporal slits and the corresponding interferences. Ex- 
perimental observations have been carried out with ul- 
tracold neutrons Q and with ultracold atoms |3, 0, 01 ■ 
Recently, it has also been observed for electrons in a dou- 
ble temporal slit experiment J\. The "shutter problem" 
has indeed important applications, as it is the modeliza- 
tion of turning on and off a beam of atoms as done e.g. in 
integrated atom-optical circuits or a planar atom waveg- 
uide ;81|, and may allow us to translate the principles 
of spatial diffractive light optics to the time domain for 
matter waves 0. Besides, it provides a time-energy un- 
certainty relation liol [T]| which has been verified exper- 
imentally for atomic waves, by realizing a Young inter- 
ferometer with temporal slits 0, |H • The Moshinsky 
shutter has also been discussed with the Wigner func- 
tion and tomographic probabilities |l2l. fo r relativistic 
equations 

B m El, with dissipatioVjil in rela- 
tion to Feynman paths |l7j . or for time dependent bar- 
riers 01, including the time dependence of states ini- 
tially confined in a box when the walls are suddenly 
removed |23| . With adequate interaction potentials 
added to the model, it has been used to study and char- 
acterize transient dynamics of tunnelling matter waves 
Ulillla^iJl, J5.,26,1I], and the transient response 
to abrupt changes of the interaction potential in semi- 
conductor structures and quantum dots . 



Within the source boundary approximation, in which 



the form of the wave is imposed for all times at a source 
point or surface, many works have been carried out 
within the field of neutron interferometry ,,3[L 31], con- 
sidering also a tr iang ular aperture functions 1331, at om- 
wave diffraction |33l |. tunnelling dynamics |l4l 1251 Is^ 
I35L I3R 37, 38], and absorbing media, in which an ul- 
trafast pe ak- p ropagation phenomenon has been recently 
described 
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The importance of pulse formation is nowadays en- 
hanced due to the possibilities to control the aperture 
function of optical-shutters in atom optics, and to the 
development of atom lasers. For such devices some of 
the first mechanisms proposed explicitly implied periodi- 
cally switching off and on the cavity mirrors, that is, the 
confining potential of the lasing mode. As an outcome, a 
pulsed atom laser is obtained. Much effort has been de- 
voted to design a continuous atom laser, whose principle 
has been demonstrated using Raman transitions js^, Eol ■ 
With this output coupling mechanism, an atom initially 
trapped suffers a transition to a nontrapped state, re- 
ceiving a momentum kick during the process due to the 
photon emission. These transitions can be mapped to 
the pulse formation, so that the "continuous" nature of 
the laser arises as a consequence of the overlap of such 
pulses. There is in summary a strong motivation for a 
thorough understanding of matter-wave pulse creation, 
even at an elementary single-particle level. 

In this paper we aim to describe the characteristics 
of one-dimensional, matter-wave pulses within the ap- 
proximation in which interatomic interactions can be ne- 
glected, as it is usually the case in standard low-pressure 
atomic beams. This will be useful as a reference and first 
step to consider the interacting case later on. While the 
most idealized pulses have been studied in several works 
since the seminal paper of Moshinsky P| , some aspects of 
realistic pulse production have not been examined with 
enough detail yet, such as the effects of statistical mixing 
in the preparation beam, or the fact that optical shutters 
do not have absolutely sharp edges, and could be taylored 
for designing the pulse features. Other main objective 
of the paper is to compare two ways to model boundary 
conditions: as a standard initial value problem where the 
wave is specified in the preparation region at time i = 0, 
or according to the "source approach" in which the wave 
function is given for all times at the source point. 
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The organization of the papers is as follows: In Sec. 
^we review the Moshinsky shutter, whose main feature 
-diffraction in time- is quantified by means of the fringe 
visibility in Sec. IIIII Section IIVI deals with the evolu- 
tion of finite pulses; In Sec. the exact results are 
compared with the usual "source" approximation, and, 
finally, in Sec. IVII the aperture function is modified for 
different types of apodization and the time-energy uncer- 
tainty product is evaluated. We also consider periodical 
shutter apertures and find analytical solutions for the 
wavef unction. 



II. MOSHINSKY SHUTTER 

In this section we review the "Moshinsky shutter" 0, 
ITol | and discuss its physical interpretation. Consider a 
plane wave impinging from the left on a totally absorbing 
infinite potential barrier located at the origin (shutter) 
which is suddenly turned off at time zero. For t = one 
has the initial condition 

V'(x,t = 0) = e*?'^/'"'e(-x), 

where Q is the step function, whereas for t > the wave- 
function can be written using the free-particle Green's 
function. 
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and 



is the "w" or "Faddeyeva" function 141|,|42|. The contour 
r_ goes from — oo to oo passing below the pole. 

The associated "density" ^ {tpix, i)p exhibits diffraction 
in time, namely, a characteristic oscillation in time, and 
also in the space domain, see Fig. Q in contrast with 
the simple step function for a spatially homogeneous en- 
semble of classical particles with momentum p released 



^ Actually |i/'(a;, t)p is a relative density with respect to the asymp- 
totic, long-time, stationary value. 




FIG. 1: Diffraction in time seen in P{0), Eq. lO, for an 
initial cut-off plane wave state. 
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FIG. 2: Cornu espiral. The points associated with the max- 
imum and minimum of the largest fringe are plotted. They 
correspond to a maximum and a minimum of the distance to 
(-1/2,-1/2). 



at time i = 0. Thanks to the relation between the w- 
function and Frcsnel integrals |4^. 
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the solution may also be rewritten as 



\^{x,t)\^ ^ p{e) 



S{9) 



C{9) 



(3) 

and therefore be mapped onto the Cornu espiral, which 
is plotted in Fig. [21 The "density" can then be read 
as half the distance from the point (—1/2, —1/2) to any 
other point of the spiral; and the origin corresponds to 
the classical particle with momentum p released at time 
t — from the shutter position. 

Suppose now that the shutter is closed suddenly 
at time r. The integrated "density" N+{t) := 
/q dx \'ip{x,T)\'^ grows continuously with r and may 
reach arbitrarily large values, obviously greater than one. 
Moreover, \'ip{x,t)\'^ is dimensionless, so N+ cannot be 
taken as a probability in any case. This poses the ques- 
tion of physically interpreting the mathematical results. 
Clearly we face analogous problems to interpret a plane 
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wave or stationary scattering states, which are not in 
Hilbert space. The solution is then found similarly, by as- 
suming that 'ijj{x, t), except for a normalization factor, is 
just a component of a normalizable state in Hilbert space, 
to be determined by the experimental preparation setup. 
A consistent interpretation for several -noninteracting- 
particles requires the appropriate quantum symmetriza- 
tion ,43], which will be treated elsewehere. In the present 
paper we shall limit ourselves to the simplest case an con- 
sider only one-particle states. 

For a shutter with a non-zero reflectivity, the cut-off 
plane wave of opposite momentum has to be taken into 
account, so the initial state takes the form 

^^^\x, t^O):^ {e'P/^'^ + 7^e-''P/'"'^)e(-a;). (4) 

The cases TZ = {1,0, —1} turn out to be the most rele- 
vant. We will refer to 7?, = — 1 as sine initial conditions, 
which arise when the shutter is totally reflective; and to 
TZ — 1 as cosine initial conditions. This last case may 
seem to be hard to implement but, surprisingly, it is re- 
lated to the usual source boundary condition. 



*e(i), tpt{x>0,t<0)^0, (5) 



where uj = p^/{2mh), and is more tractable mathemat- 
ically than the sine case. As we shall see, ^pp{x,t) = 



III. VISIBILITY OF MAIN DIFFRACTION 
FRINGE 

Let us examine flrst the time evolution from a cut-off 
plane wave initial condition, i.e., 7?. = 0, as in Eqs. 
and (jSJ. The motion of points of constant "density" is 
obtained by imposing a constant value of 9, according 
to Eq. ®. In particular, for the maximum and nearby 
minimum, see Fig. |21 



Xmax(t^ pt/m \ 9jji(ix: 

" m 



irfit 



Xniinii^ — pt I Tfl w ^min^ 

V m 

where Omax — 1.217 and Omin = 1.872 are universal for 
TZ = 0. Since the probability is exclusively ^-dependent, 
see Eq. Pmax — 1-370, and Pmin = 0.778, indepen- 
dently of time, mass or momentum, and therefore the 
fringe visibility, defined as 



V(T) = 



Pi. 



Pl^^ max -^1^ 



is concluded to be universal (equal to 0.276 for 7^ = 0), 
and of course time- independent. 

Other important fringe feature is the width of the main 
peak. It can be computed from the intersection with 
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FIG. 3: Visibility as a function of temperature and the reflec- 
tivity of the shutter for a beam of Argon with the distribution 
in Eq. Q. pc/m = 10 cm/s, t = 20 /is. 



the classical probability density (for a cut-off beam of 
particles released at time with momentum p) P, Il0j| , 



Ax ~ 0.85 
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Next, we consider an initial state given by a statistical 
mixture, 

p{t = 0) :^ / dp/(p)|^f )(t = 0))(^f Hi = 0)1, (6) 



and focus our attention on the effect of temperature and 
evolution time on the diffraction pattern. Let us as- 
sume, by integrating out the y,z-degrees of freedom of 
a Maxwell-Boltzmann momentum distribution j44j . that 
the Gaussian momentum probability density 



fMsip) 



1 



{p — pc 



(7) 



holds for the momentum in x direction. The introduction 
of a non-zero mean average momentum p^ implies that 
the atomic ensemble, say in a magneto-optical trap, is 
launched along the a;-axis with a velocity Pc/nn. This is 
realized, for example, in the "moving molasses technique" 
[4^ within an atomic fountain clock. We shall generally 
assume that the contribution of negative p in Eq. ([TJ is 
neghgible. 

Figure 13 shows the effect of the temperature on the 
visibility of the main diffraction fringe for statistical mix- 
tures weighted by the distribution of Eq. Q, with 
TZ = {0, ±1}. It is clear that the diffraction peaks are 
more visible for sine conditions than for the isolated cut- 
off plane wave and cosine conditions, the difference be- 
ing more noticeable at small temperatures. As expected, 
there is a suppression of the diffraction pattern when the 
temperature increases. A possible criterion for this sup- 
pression is a small value of the ratio, a, between the width 
of the fringe Ax and the separation Axc; of the classical 
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FIG. 4: Visibility as a function of the a-parameter for an 
Argon atom with a velocity Pc/m — 10 cm/s. Each curve 
is computed by varying the temperature for a given opening 
time. TZ = —1. 



arrival points corresponding to fast and slow components 
of the statistical mixture, namely, 



Ax 
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For the Gaussian distribution, Eq. |7I), one can choose 
Ap as the full width at half maximum (FWHM), that is, 
2y/2ln2mkBT, so that 



0.36 



KbTI 



According to Fig. 01 the visibility grows with a for 
small values and saturates for large a, i.e. for small times 
and/or temperatures. 

Moreover, if instead of Eq. Q, the momentum distri- 
bution of an effusive atomic beam is used HI, Hal , 



fbeam (p) 



P 



2{mkBTY 



^e(p). 



(8) 



then {p) — y^9TrmkBT/8, and the upshot is a complete 
supression of the diffraction in time phenomena for all 
times and temperatures. This fact points out the relevant 
role played by the momentum distribution, that is, by the 
experimental preparation setup being considered. 

Summing it up, the fringe visibility, at variance with 
the pure state case, becomes time-dependent for mixed 
states. The characteristic oscillations of diffraction in 
time tend to be washed out with the observation time af- 
ter opening the shutter, and also by increasing the tem- 
perature, or for the momentum distribution of an effusive 
beam. 



IV. EVOLUTION OF A FINITE PULSE 

In this section we focus on the formation and evolution 
of a matter-wave, single-particle pulse. Suppose that the 



shutter is opened at the instant t = and closed at the 
instant r. We will refer to the duration r as the "opening 
time". It is, in other words, the width of the the time 
slit. The notation V'po.HO will denote the state formed 
in that manner, at time t > t, from the initial state of 
Eq. Q with main momentum po-'^ 

For an initially pure state the time evolution of the 
pulse can be written in terms of the inner products with 
the eigenstates of the free hamiltonian which vanish at 
the origin, 



{x\(l>L 



l2mh 
irp 



sin(px/?i), E — p'^/{2m). 



They obey orthonormality and completeness relations, 



Jo 



dE\ 



As shown by Moshinsky [l^, the overlap between these 
eigenstates and the state that results from a state formed 
after an opening time r with 7?. = is 



mpfi 



w[-uo(po,t)] w[-uo{p,t)] w[-uo{-p,t)] 



pI 



2p{p-po) 2p{p + po) 



(9) 



where 



uoip,T) 



to be compared with Eq. (01. {uq{p,t) is a particular 
case of u{p,t) in which x — 0.) The overlaps for the 
cases TZ = ±1 are obtained easily from Eq. ©. 

Writing the final time of observation as < = t + At, the 
evolved pulse can then be calculated as 

/•OO 

Jo 

(10) 

This expression is not well-behaved numerically, but such 
a problem can be swiftly solved by modifying the contour 
of integration in the complex p-plane, since several terms 
can be written as w-functions (see the Appendices A and 



^ Another way to form a pulse is to remove at time t = the 
potential walls that confine a state which is initially stationary. 
For an infinite-wall square box the result is the substraction of 
two semi-infinite 1Z = —1 pulses, i.e., at least for free motion, a 
simple combination of to-functions. For applications of this i dea 
in free -motion cases or scattering configurations see llfll l2(l 123 . 

Ellll. 
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B for details), and for the rest we gain an integrand with 
Gaussian decay. ExpUcitly, for 7?, = —1, 

I 1 1 e 2th 

%o!T{x,t) = -— {[w[-u{po,t)]+w[-u{-po,m 



e *2,„n - w[uo(po,t)] ™xi 

g 2A.th 



X {w[-u{po,M)]+w[-u{-po,l^t)]} 

2tt Jr , 



w[uo{~p,t)] w[uq(p,t)] 



P + Po 



P-Po 



(11) 



where the contour passes above the poles. This result is 
so far exact. The integral term, 3^, can be rewritten by 
integrating in the u variable, 



y 



le 2AtR 

2^ 



due-"" 



u-u{po,t) u-u{-po,t) 
see the Appendix A, with 

1 + i It mx 



(12) 



mh At 

A way of dealing with this kind of integral to excellent 
accuracy was discussed by Brouard and one the authors 
in [231 . It consists in extracting the singularities of the 
integrand, which we shall denote as g{u), in the following 
way 



h{u). 



Ai is the residue of g{u) at the pole Ui and h{u) is an 
entire function which admits a power series expansion 
convergent in the whole u-plane, 



h{u)e " du 

) 

fe(u = o) + y ~\^,' fe(^")(z. = o) 

n=l ^ ' 



Retaining just the first term of the previous expansion 
leads to 



In particular. 



'ITTW^ — Ui) H 



+ v^.9(u = 0). 



y 



g 2At/i 



{w[uo(po,T)]w[-u{po,At)] 



+ w[-uo(po,t)]w[-7x(-po. At)]} 



20F t u{pa,At) 
w(-w)-w[-uo(po,t)] 



(13) 



u{~PQ,At) 

which we shall use later on in Sec. to compare the 
initial value problem and source approaches. 



A. Purity of the chopped state: coherence creation 
from incoherent mixtures by using small opening 
times 

It is possible to "create coherence" from statistical mix- 
tures, simply by means of small opening times. This was 
seen experimentally with cold atoms by Dalibard and 
coworkers @. To quantify this effect we use the purity, 
defined by 

Vt :=Trp2,^,(t), 

where pN.rit) is the normalized density operator at time 
t for an opening time r. Our starting point is the unnor- 
malized state 



Prit) := /dp/(p)|4")(t))(^(")(t)|. 



(14) 



Thanks to the invariance under cyclic permutation of the 
trace and unitarity of the time evolution, the purity is 
invariant for t > t, so it can be calculated in particular 
at t — T, 



1 



^2 / / dpdp7(p)/(p')l(V'i,°)l<,V)l^ 

where N is given by 

N = TTpr = I dp/(p)(^W|^W). 

Considering the case of an effusive atomic beam, /(p) = 
fbeamip), Eq. © , Fig. shows that the purity Vt de- 
creases as a function of the opening time and temper- 
ature. For short opening times a highly coherent state 
is created in spite of the mixed preparation state. This 
result, in the limit t ^ 0, may be understood as follows: 
As long as a (pure state) wavefunction vanishes in x > 
at time t = Q, the wavefunction for t > t' > can be 
written as 46] 

ipix^t) = [ dt'K+{t,x:t' ,x' ^Q)il}(x' = 0,t'), 
Jo 

where 

K+{t,x;t',x' = 0) ^ —-^go{x,t;0,t') 
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FIG. 5: Variation of purity versus the opening time for dif- 
ferent temperatures. An effusive beam of Argon atoms is 
considered, see Eq. ©. TZ = 0. 



By performing a Taylor series expansion in t' around t' — 
0, one finds that 

t) = tK+{t, x; 0, 0) + 0{t^) 

for ^{x = 0,i' = 0) 1. Since ^l!^,l{x = 0,t' = 0) = 1 
for all p, then to first order all the momenta that form 
the preparation mixture give the same contribution, in 
other words, |'0p,r(i)) in Eq. p4|l becomes independent 
of p, and thus the resulting state tends to be indepen- 
dent of the initial distribution f{p), and totally coherent 
(pure) in the limit of small opening times. For the case 
of the Gaussian distribution in Eq. Q, the same qual- 
itative behaviour (monotonic decay from 1 in the same 
time scale) is observed. 



V. COMPARISON BETWEEN "SOURCE" AND 
"INITIAL VALUE PROBLEM" APPROACHES 

Golub and Gahler [s^ tackled the diffraction in time 
problem in two dimensions using the Green's function 
formalism, working within the Fraunhoffer limit, and 
imposed source boundary conditions. This approach 
was later thoroughly implemented for different combina- 
tions of time and space rectangular slits by Brukner and 
Zeilinger fs^l. Other works have used the simple source 
boundary conditions of Eq. O in one dimension, see 
[38j |. and references therein. In |46j| the relations between 
source boundary conditions and the initial value problem 
(IVP) were examined and equivalence conditions estab- 
lished. 

In order to compare pulses formed by source or IVP ap- 
proaches let us introduce the window function X[o,t](0 = 
Q{t)0{T — t) and impose the source boundary condition 



V';„^.(a; = 0,0 = e-*"»*X[o,r]W, 



(15) 



all t < 0) with Fourier transform 

^;^„,(x = 0, c.) = ^ \1 - o)^ 



At a later time one finds, for a; > 0, 



duj 



pipx /h—iut 



' ' 27r U -Uq 



where Imp > 0. This can be written in the complex p- 
plane by introducing the contour of integration r_|_ which 
goes from — oo to oo passing above the poles. 



KA^,t) . - j^^ dp2p-^-^[l - eH^-)^ 

= [ dp ( ^_ + ^_ ) e*P"^/''-*5^fl-e^("""'')^ 
27rJr+ \P + Po P-PoJ 

The integral can be carried out exactly using Appendix 
A, 



e 2tfi 



{w[~u{pq, t)] + w[-u{-pQ, t)]} 



{w[-u{po, At)]+w[^u{-~po, At)]}, 



(16) 

so the time evolution is straightforward and given by a 
difference of w- functions (in fact, expressing the win- 
dow function as a difference of step functions X[o,r] = 
Q{t) — Q{t — t), the solution can be foreseen to be given 
by the difference of one source open at t = and a sec- 
ond one at a later time t = t, with the associated phase 
taken into account). The nice thing about this approach 
is that it allows a clear and intuitive interpretation of the 
different contributions to the total wavefunction, and it 
is easy to implement complicated combinations of time 
and space windows. The flip side is that deviations from 
the exact result of the initial value problem arise, spe- 
cially for the back part of the wavefunction cut by the 
chopper, except in the case TZ = 1. Although the calcu- 
lation is somewhat tedious -see Appendix C-, it can be 
proved that for cosine initial conditions, the wave func- 
tion that results from a shutter opened between and r, 

''Pp^jlrix, t), is exactly equivalent to using source boundary 
conditions in that interval, i.e., to Eq. I|16() . The differ- 
ence with the solution for the reflecting shutter {TZ = —1) 
can be written as 



4~'Ax,t)-i;<^2ri^,t) 



"pa.T y-^T-/ -rpQ.ry-J-/— 2 

X {w[uo{po, t)] - w[-~uq{po, t)]} w[-u{-po, At)] 
ie-i'^tn I w{v) — w[uo{po,t)] 



(the wave function is also assumed to vanish at a; > for 



V uipo,At) 
w{-v)-w[-uo{pq,t)] 
u{-po,At) 



(17) 
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FIG. 6: Overlap probability between the wavefunction of one 
Argon atom associated with sine and cosine initial conditions 
(continuous line), cut-off plane wave and sine conditions (cir- 
cles), and cut-off plane wave and cosine conditions (dashed 
line), as a function of the phase S/h :— rp^/ {2mh), see Eq. 
itTHl . 



where we have used Eqs. and (|16|l . The overlap 

probabihty between pulses generated for different values 
of TZ and fixed r and po , 



)\|2 



(18) 



is independent of the evolution time and can be cal- 
culated exactly at t = r. It turns out to be a function 
of the adimensional phase S/Ti := rpo/(2m?i), see Fig. 
El The overlaps tend to one as t ^ 0, oo. Remarkably, 
there is a lower bound which is 93.56% between sine and 
cosine— source conditions. We may conclude that in gen- 
eral the result of using source boundary conditions is not 
at all a bad approximation to the solution of the initial 
value problem with reflecting or absorbing conditions, 
except when very accurate results are required. 



VI. APODIZATION OF MATTER- WAVE 
PULSES 

Apodization is a technique long used in light optics to 
avoid the diffraction effect in several devices. Simply put, 
it consists in suppressing high frequency components by 
smoothing the aperture function (or "window" for short) 
of the pulse [43 • The price to pay is a slight broad- 
ening of the energy distribution. In this section such a 
technique is applied to matter wave pulses within the 
source approach. In particular, the case of a sine type 
aperture function can be analytically solved in terms of 
w-functions. Let us consider in Eq. (|15|l . instead of the 
sudden, rectangular shape, a sine aperture function. 



x[o,r]W -sin(m)e(i)e(T-t) 



(19) 



with r = ttJI ^ . (Note that the "sine" aperture function 
has nothing to do with sine initial conditions.) Then, 



defining u:± = ujq ±0. and the corresponding momenta 
p± = ^2mfiu}±, with the branch cut taken along the 
negative imaginary axis of uj±, it is found that 



1 



UJ ~ UJ- 



and the wavefunction is (we shall distinguish apodized 
pulses from the rectangular aperture pulses used so far 
by means of a tilde) 



-E 

27r ^ 



a 



a=± 



duj e-*'^*+»P^/'i 



1 — ei{'^-'^€,)T 



—— a dpe 



a = ± 
1 



ipx /h I ^ — iujt ^—iujAt — iuj, 



1 



P+Pa P-Pa, 



- 2^ a e— {w[-u{po,,t)]+w[-u{-p^, t)]} 



a=± \ 

_ e '2™R+2AtK{i„[_u(p^^ At)]+w[-u(-pa, Ai)]} 

As it is clear from the last line, the apodization with the 
sine aperture function is equivalent to the introduction 
of two different momenta components whose values are 
related to the length of the pulse. Pictorically, a sine win- 
dow pulse is equivalent to the coherent difference of two 
simultaneous pulses produced with sudden rectangular 
windows from speeded up and decelerated sources. 

Similarly, other window functions, as the ones shown 
in Fig. 13 can be worked out. For instance, the Hanning 
aperture function, 

x[o,r]{t) = sin\nt)e{t)eiT-t) 

= i(l-cos2m)e(i)e(T-i), (20) 



after introducing pp = y^2mh{uJo ± 2fl) with /3 = ±, 
leads to 



V'p'o,r(a;,t) = 



1 



/3=± 



where the effect of the apodization is to subtract to the 
pulse with the source momentum two other matter-wave 
trains associated with pp, all of them formed with rect- 
angular aperture functions. 

A related configuration wich can be solved analytically 
in terms of ly-functions is the source apodized periodi- 
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FIG. 7: Different aperture function studied. The smoothness 
at the edges increases in this order: rectangular (continu- 
ous Une), sine (dotted Hue), Hanning (dashed) and Blackman 
(long dashed hne). 
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FIG. 8: Front part of the probabihty density for the peri- 
odic Hanning function a.t t = 1.2 ms, with an aperture time 
r — lOfis. The circle, diamond and square mark respectively 
the classical position of the slow (p-), original (po) and ac- 
celerated (p+) components, po/m = 10 cm/s. 
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FIG. 9: Probability density of a pulse of ultracold Argon 
atoms {po/m — 10 cm/s), with r = 10 /is registered at t = 
200 fis for different window functions: rectangular (continuous 
line), sine (dotted line), Hanning (dashed line) and Blackman 
(dot-dashed line). 
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FIG. 10: Detail of the sidelobes of the previous figure. 



cally to create successive pulses. This can be done ex- 
tending Eq. (|20|l to all positive times. The result is 

~ e 2th 

V'po,r(2^.*) = -^—{M~'uiPO,t)] +w[-u{~po,t)]} 

1 

%=± 

Note again the intervention of three different momenta. 
A consequence is that the w-function for the fastest com- 
ponent will eventually separate spatially as a smoother 
-rather than pulsed- forerunner, see Fig. |S1 

Also suitable for suppression of the sidelobes is the so- 
called Blackman aperture function, 

X[o,t] (t) = 0.42 - 0.5 cos(2rJi) + 0.08 cos(4rJt) . 

Using the defined above and the momenta 

= \/2mh{LUQ ± 4r2), 



with f3,"f ~ ±, one can write the time evolved wavefun- 
tion in terms of twenty w- functions, 

V^;„,,(x,t) = 0.42^;„^,(a;,i) - 0.25 ^ 

/3=± 

+ o.o4^^;^,,(x,t), 

SO, once again, the smoothed shutter is related to a com- 
bination of different pulses from rectangular time slits. 

The probability density for a single pulse and the de- 
tails of the secondary diffraction peaks for four differ- 
ent aperture functions arc illustated in Figs. 1^ and [TUl 
where all states are normalized to one. As expected, 
the smoother the time window function, the higher the 
suppression of the sidelobes. Note that the Blackman 
apodized pulse is slightly more advanced that the other 
ones due to the faster components induced by this par- 
ticular apodization. 
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FIG. 11: Uncertainty product, AEr, for a pulse of ultracold 
Argon atoms (po/m = 10 cm/s); sine aperture function and 
Ai? = FWHM (solid line) ; sine aperture function and A_E = 
cr_B (dotted line); rectangular aperture function with AE = 
FWHM (dashed line). 



A. Time-Energy uncertainty relation 

A time energy uncertainty relation was discussed for 
the shutter problem by Moshinsky, who calculated the 
energy distribution of a particle after closing the shutter 
[Tol| . Using cosine conditions with energy E and taking 
the overlap of the wavefunction with the free particle 
eigenstates which vanish at the origin, he calculated the 
energy distribution when the shutter was open a time r, 



P{E,E',t)=N^/E- 



-sin\E - E')T/2h] 
{E - E'Y 



A/" being the normalization constant. Then it was con- 
cluded that, for some measure of the energy width ^ 
(see [T]| for a comprehensive review of time-energy un- 
certainty relation), 

i^E T ~ n. 

In fact, if one introduces the root of the variance ue 
as the measure of energy spreading, Ai? = aE^ such an 
uncertainty relation cannot be established since ue di- 
verges. This is in consonace with the non-existence of 
the average displacement, (x), for sharply cut waves [48j| . 
One can calculate also the energy distribution for the case 
of a sine-aperture function and source boundary condi- 
tions, Eq. H19|) . analytically. It turns out to be slightly 
broadened, but the variance is now well-defined as a con- 
sequence of the supression of high frequency sidelobes. 

Alternatively, as done in 6] , the FWHM can be taken 
as the measure of energy spread. One observes that in 
both cases, sudden and smooth shutter, both measures 
of the spread decrease with the opening time. The uncer- 
tainty product is represented versus the opening time in 
Fig. ^] Note the two regimes, first a linear growth and 
then saturation, separated by the characteristic time re- 
quired for a classical particle to travel a de Broglie wave- 



T 



hm 



VII. CONCLUDING REMARKS 

Starting with a brief review of the Moshinsky shut- 
ter problem, we have examined several one-dimensional, 
matter-wave pulse characteristics such as the visibility of 
the diffraction-in-time fringes versus evolution time and 
temperature, or the creation of coherence by small aper- 
ture times. Moreover, a long standing missed comparison 
among several initial conditions used in previous works, 
which differed on the reflected components in the prepa- 
ration region, has been carried out. We have studied sev- 
eral apodizations of the aperture function, and obtained 
analytical expressions for the time evolution of the corre- 
sponding pulses, also in the periodic case. The effects of 
apodization to suppress secondary diffraction peaks and 
in the energy width have been discussed. 

The present work is the first step towards a more am- 
bitious objective, namely, tayloring the pulses for partic- 
ular needs, including the case in which interatomic in- 
teraction is important. Applications for laser atoms, in 
particular, will require in general to consider non-linear 
effects. Physically, optical shutters formed by effective 
potentials due to detuned lasers offer a unique oportu- 
nity, not available with mechanical shutters, to control 
the formation and subsequent behavior of the pulse. 
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APPENDIX A: INTEGRALS 

Consider the following integral along a contour which 
goes from — cx) to oo passing above the pole at pq, 



dp- 



-iap -\-ibp 



P-Po 



-, a > 0. 



The saddle point is at ps — b/2a. By completing the 
square, one is lead to introduce the variable 



u = u{p) — —j=- ( Vap — 



V2 



Deforming the contour to go along the steepest descent 
path from the saddle, which in the u-plane is at the ori- 
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I = e*4a / du 



r„ u- u{po) ' 

If Co denotes a counterclockwise circle around the pole 
m(po)j then r„ = (—00,00) U Cq if the pole has been 
crossed {Im{u) > 0) and the real line of u otherwise. 
We can inmediatly recognize this integral as a Faddeyeva 
function since the following equalities hold thanks to 
Cauchy's theorem, 

— / du 



in u - u{po) 

du- 



2g-["(po)]^ 



in Jy_ u — u{pq) 

(po)] - 2e-["(f«)l' - -w[-uipo)], 



w\u 



where r_|_ and r_ go from —00 to 00 passing above and 
below the pole, respectively. One concludes that 

X = —iire''^'' w[—u{po)]. 



APPENDIX B: SINE INITIAL CONDITIONS 



where the relation 

2 

w{z) + w{—z) = 2e^^ 
can be applied leading to 

/OO ,2 



(B2) 



X <; e 's^fi - i(;[uq(p, t)] - e *2™k 



with 



p' + p p' — p 



/•LXJ ,0 



(B3) 



P - P 



p' +p 



The resulting integrals, but for the term are of the 
form solved in Appendix A. The final result is Eq. (|ll|l . 



After opening the shutter and closing it at time r, the 
wavefunction for a reflecting shutter, TZ = —1, is given 

by 

V'(-i) {x, t) = V^o) (x, t) - i>%{x, t), 

which is a difference of w-functions. 
Since, using Eqs. I© and (fTU)) . 

tpi°l{x, i) = — / dp'p' e-'"^ sm{p'x/h) 
27r J^^ 

\w[-uo{p,t)] w[-~uo{p',t)] w[-uo{-p',t)]\ 



p'2 _ p2 2p'{p' - p) 2p'{p' + p) 



(Bl) 



the wavefunction can be written as 

dp'e-^'^ie'P'''/^ - e~'v^in-^ 

-00 

^ j w[Mo(y,r)]-w[~uo(y,T)]-w[-Mo(p,T)]+w[Mo(p, r)] 
\ p' +p 

^ w\-Uq{p' , t)\-w\uq{p' , t)]-w[-mo(p, r)]+w[Mo(p, r)] [ 

P' -P J 

Using the symmetry under p' —p' the previous expres- 
sion is simplified to 

/OO ,2 
-00 

^ I wjuoip' ,t)]-w[-~uo{p' ,t)]~w[~uo{p,t)]+'w[uo{p,t)] 
I p' +p 



APPENDIX C: COSINE INITIAL CONDITIONS 

The solution to the Moshisnky shutter with cosine ini- 
tial conditions, TZ = 1, 



(a;, t) = {^«[-^(P, i)] + w[~u{~p, t)]}, 

is immediately computed from Eq. (Q. The same result 
is obtained using the source boundary conditions of Eq. 
(O and following the steps leading to Eq. H16|) . 

In fact the same agreement can be extended to finite 
opening times r. The procedure to obtain ippHriXji) is 
very much the same as for the sine case. Using 

^^]lix,t)=^/jl,"}{x,t)+i:^°l,{x,t), 

the expressions for the ip'fp ^ (t) in Eq. ljBl|l , the identity 
in Eq. HB2|) . and t = t + At, one finds 



4]l{x,t) = ^ / dp' [e- 



/h ~ip' X /h 



p' + p p' — p 



+ 



w[-Uo{p' ,t)]-w[uo{p' ,t)]-w[-uo(p, t)]+w[uo{p, t)] 



P -P 



Due to the symmetry under p' — > —p' the term with 
_g-ip x/h gq^g^j ^YiQ with e^P and the result- 
ing integral can be carried out by completing the square 
and deforming the contour of integration in the complex 
p-plane as described in Appendix A, so new Faddeyeva 
functions can be identified. This leads to the same ex- 
, pression found in Eq. pt)|l . 
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